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Abstract 

We present a study of the critical phenomena around the quantum critical 
point in heavy-fermion systems. In the framework of the S = 1/2 Kondo 
lattice model, we introduce an extended decoupling scheme of the Kondo 
interaction which allows one to treat the spin fluctuations and the Kondo 
effect on an equal footing. The calculations, developed in a self-consistent one- 
loop approximation, lead to the formation of a damped collective mode with 
a dynamic exponent z = 2 in the case of an antiferromagnetic instability. The 
system displays a quantum-classical crossover at finite temperature depending 
how the energy of the mode, on the scale of the magnetic correlation length, 
compares to ksT. The low temperature behavior, in the different regimes 
separated by the crossover temperatures, is then discussed for both 2- and 
3-dimensional systems. 

PACS numbers: 75.30.Mb, 71.27.-ha, 75.50.Ee, 71.28.-hd 
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I. INTRODUCTION 



A central issue in the thermal properties of heavy- fermion compounds is the discovery 
of a non-Fermi liquid behavior in systems close to an antiferromagnetic quantum phase 
transition |T| . Such a behavior has been discovered in a series of compounds containing 
cerium or uranium, for example, CcCuq-xRx {R=Au, Ag) 0, Celris, CePd2Si2 0, 
CeNi2Ge2 0], Ui-xY^Pds ^ and CexLai_xRu2Si2 0- In CeCu^,gAuo,i, the specific 
heat C depends on T as C/T ~ —ln{T/To), the magnetic susceptibility as x ~ l — aVT, 
and the T-dependent part of the resistivity as Ap ~ T (instead of C/T ~ x ~ Const 
and Ap ~ in the Fermi liquid state). Most interestingly, the breakdown of the 
Fermi liquid behavior can be tuned by alloying (chemical pressure) or by applying a 
hydrostatic pressure or a magnetic field. The origin of this NFL behavior is a problem 
of current debate. 

Different scenarios have been proposed in order to elucidate the non-Fermi liquid 
(NFL) behavior at the quantum critical point (QCP). On the one hand, there are sce- 
narios involving strong disorder to account for the non-Fermi liquid properties, including 
(i) a distribution of the Kondo temperature P{Tk) induced by disorder following a 
distribution of the local density of states, and (ii) the possibility of the formation of a 
Griffiths phase at the quantum critical point On the other hand, starting from a 

clean or rather weakly disordered system, there is a scenario (iii) refering to the prox- 
imity of a quantum critical point in a theory of itinerant antiferromagnetism |T0|-p]7| , in 
which the Fermi surface is coupled to critical collective antiferromagnetic excitations. 
An alternative possibility has been recently put forward |]TB[ on the basis of an analysis 
of the neutron scattering measurements in CeCu^,gAuo,i. The unusual spin dynamics 
that has been observed suggests (iv) the existence of critically screened local magnetic 
excitations due to the formation of local moments. Finally, it has been proposed (v) a 
variant of the quadropular Kondo effect [|l9l in which the interactions of the rare earth 
(actinide) internal degrees of freedom with the conduction electrons are reduced to a 
single impurity multichannel Kondo effect. As it stands, the situation remains highly 
controversial. 

In this paper, we adopt the point of view (iii) of the proximity of a quantum critical 
point. The effect of a non zero temperature on quantum critical points is a long-standing 
problem encountered in itinerant magnetism. We refer to the papers of Hertz and 



Millis using renormalization group techniques, and of Moriya introducing a 
self-consistent renormalization (SCR) theory of spin-fluctuations. These theories yield 
a rich phase diagram with, in addition to the ordering temperature T^, a series of 
crossover temperatures separating different regimes of behaviors. They are well adapted 
for the treatment of spin fluctuations as present in the single-band Hubbard model but 
essentially do not account for the Kondo effect. 

Indeed, the model which is believed to describe the physics of heavy fermions is the 
Kondo lattice model (KLM) [E^ derived from the periodic Anderson model (PAM) in a 



given limit. The large N expansions ||23| which have been carried out for these models 
(where N represents simultaneously the degeneracy of the conduction electrons and of 
the spin channels) give a good description of the Kondo effect. But, unlike the above- 
mentioned spin-fluctuation theories, they fail to account for the spin-fluctuations since 
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the RKKY interactions only appear at the order ||24|| . 

Hence our motivation to set up an approach of the S = 1/2 Kondo lattice model 
{N = 2) that enlarges the standard expansion theories up on the spin-fluctuation 
effects. Our previous paper pSf , presented such an approach by performing a generalized 
Hubbard-Stratonovich transformation on the Kondo interaction term of the Hamiltonian 
which makes three types of fields (magnetization densities and Kondo-type) appear on 
an equal footing. The dynamical spin susceptibility was derived in a one-loop expansion 
associated with the gaussian fluctuations of the magnetization density fields around 
their saddle-point values. 

With the aim to describe the critical phenomena around the quantum critical point 
of the heavy-fermion systems, we propose to push the treatment of the S = 1/2 Kondo 
lattice model in a self-consistent one-loop approximation in complete analogy with the 
SCR theory of spin fluctuations developed for the Hubbard model pTip6| , p7| . As in itin- 
erant magnetism, we will show how the system displays a quantum-classical crossover 
at finite temperature depending whether the temperature is lower or higher than the 
characteristic energy-scale of the collective mode associated with the magnetic insta- 
bility. The low temperature behavior in the different regimes essentially depends on 
the value of the dynamic exponent z of the mode and on the dimensionality d of the 
problem. In the case of an antiferromagnetic instability, z is found to be equal to 2. 
Heavy-fermion systems are usually believed to be three-dimensional. However, a recent 
proposal ||28| based on the neutron scattering data obtained in CcCuq-xAux stipulates 
that the critical magnetic fluctuations which the quasiparticles are coupled to, are simply 
two-dimensional. So both situations d = 2 and d = 3 will be considered. 

The rest of the paper is organized as follows. In part II, we introduce the self- 
consistent one-loop approximation method for the Kondo lattice model using a gen- 
eralized Hubbard-Stratonovich transformation on the Kondo interaction term. The 
approach combines both aspects of the 1/N expansion and of the self-consistent renor- 
malization theory of spin fluctuations. The calculations, presented in the functional 
integral formalism, closely follow the presentation made by Hertz and Klenin for 
the Hubbard model. In part HI, we extract the general expression of the dynamical 
spin susceptibility extending the results of the previous paper |2^ to a self-consistent 
treatment. Part IV gives a discussion on the nature of the damped collective mode asso- 
ciated to the proximity of the magnetic transition with a dynamic exponent z = 2'm the 
antiferromagnetic case. In part V, we derive the resulting quantum-classical crossover. 
This occurs at a finite temperature T, which depends on how the energy of the mode, on 
the scale of the magnetic correlation length ^ compares to fc^T. The low temperature 
behavior of the system in the different regimes, separated by the crossover temperatures, 
is then discussed in part VI. Concluding remarks follow in section VII. 



II. SELF-CONSISTENT ONE-LOOP APPROXIMATION TO THE KONDO 

LATTICE MODEL 

We consider the Kondo lattice model (KLM) as given by a periodic array of Kondo 
impurities with an average number of conduction electrons per site < 1. In the grand 
canonical ensemble, the Hamiltonian is written as 
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ka 



^ ka 



(1) 



in which S/j represents the spin (5 = 1/2) of the impurities distributed on the sites 
(in number Ns) of a periodic lattice; c|.^ is the creation operator of the conduction 
electron of momentum k and spin quantum number a characterized by the energy = 

— J2 tij exp (ik.Rjj) and the chemical potential /x; r are the Pauli matrices (r^, r^, r^) 

<hj> 

and T° the unit matrix; J is the antiferromagnetic Kondo interaction (J > 0). 

We use the Abrikosov pseudofermion representation of the spin Sfi = J2 fL'^'aa' fia' ■ 

era' 

The projection into the physical subspace is implemented by a local constraint 

Q. = E/i^/--i = o. (2) 

a 

A Lagrange multiplier Aj is introduced to enforce the local constraint Qi. Since [Qi, H] = 
0, Xi is time-independent. In this representation, the partition function of the KLM can 
be expressed as a functional integral over the coherent states of the fermion fields 



Z = J Vci^Vfi^dXiexp - C{T)dT 

where the Lagrangian C{t) is given by 

£(r) = Coir) + Hoir) + Hj{r) 

^Oi'^) = J2(^t^rCia + fldrfia 



(3) 



Ho{r) = J2 (^kcL^ka - f^Ns [ <^kaCka " ) + ^ XiQi 

ka V^^S ka ) i 



Following ||25[, we perform a Hubbard-Stratonovich transformation on the Kondo 
interaction term Hj{t) which makes the fields $* and ^j-., appear. We get 



with Sc, = E c\„T„„,Ci„i. 



(4) 



with 
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C'{t) = Co{t) + H,{t) + H'j 



H'Ar) = j:(cl ft 



'^O ^ I ' cru' 




with Js = J/^ and Jc = ^/3. 

The saddle-point solution is obtained when neglecting the space and time dependence 
of the different fields. Details of the solution are given in Appendix A. It gives rise to the 
band structure illustrated in Figure 1 with the formation of a Kondo or Abrikosov-Suhl 
resonance pinned at the Fermi level and split by a hybridization gap defining two bands 
a and [3. The saddle-point solution transposes to N=2 the large-N results obtained 
within the slave-boson mean-field theories (SBMFT). 

Then we consider the self-consistent one-loop approximation beyond the saddle-point 
solution in the magnetically-disordered phase ^/g = ~ 0- Following Read and Newns 
we take advantage of the local U(l) gauge transformation of the lagrangian C'{t)\ 
rjexp(i6'j); f-exp{i9i); Aj —>■ + i dOi/dr. We use the radial gauge in 
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which the modulus of both fields and $* are the radial field rj, and their phase 9i 
(via its time derivative) is incorporated into the Lagrange multiplier Aj which turns out 
to be a field. Use of the radial instead of the cartesian gauge bypasses the familiar 
complications of infrared divergences associated with unphysical Goldstone bosons. We 
let the fields fluctuate away from their saddle-point values : rj = tq -|- 5rj, A, = Aq + ^Aj, 
^f. = 5^f. and = ^^c- ■ After integrating out the Grassmann variables in the partition 
function in Equation @, we express Z as a functional integral over the bosonic fields 



in, Ai 



where the effective action is 

+ nJcY.^l + JsY. ih ■ ^c. + Nsifin, - Ao)], 

i i 

with : 



(5) 



(6) 



(CTO + JcSri)Saa'Sij 



(7) 



As in the conventional approaches to critical phenomena, we have integrated out the 
fermions and thereby reduced the problem to the study of an effective bosonic theory 
describing fluctuations of the ordering fields. But contrary to Hertz ||2^ and Millis 



11 1 , who considered a single field corresponding to the magnetization density of one 



type of electron, our approach, starting from a different microscopic model, extends the 
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description to the case of several fields i.e. the Kondo fields (n, Aj) on the one hand, 
and the /- and c- magnetization densities on the other. 

We then propose to solve this field theory in a variational way. We replace the full 
expression Seff{ri,Xi,^j.,^^J — S^^^j: (where Sf^j and F^°^ are the zeroth order effec- 

tive action and the related free energy) by a trial quadratic form ASyf{r, A, ^j, — 

SefM^ ^i'^/i'^ci) - S^eff in the fields Svi , 6\^, 6^f^ and 6^^. 



+ 



Si 



c—q, — v 



f-Q,-" 
c—q, — v 



+ ( SiJ,,^ K^,. ) Di-'(q, i^^) ( )]• (8) 

where the indices g, z/ stand for (q, idj^^). The elements of the matrices D^^(q, io;,/), 
d| ""^(q, icjj,) and D^~^(q, icijjy) arc determined variationally by minimizing the free en- 
ergy or, more precisely, the upper bound to the free energy according to the Feynman 
variational principle F < F^fj + (%/(r. A, - Sg^_^(r, A, ^j) where F^j^ 

represents the free energy related to -S'e//!^' ^c) ^^d (A) means the average of the 

— (2) 

quantity A{ri, Aj, ^j., over the distribution of fields exp[—ASlfj:{r, A, ^j, ^^)] i.e. 

_ (/Pr,m,Pe/,PC.-4(n,A,,e/„Cjexp[-Ag^)/r,A,^^,0]) 
(/ VnVXmfPic,) exp[-A^g)/r, A, C/, 4)]) 

After minimization, one can show that \DQ^{ci,iijjy)]rx is equal 
to {1/2)8'^ Seff{ri,Xi,^f.,^^.)/{drq^„dX-q^-^) averaged over the mentioned distribution 
of fields (and equivalent expressions for the other elements). Explicitly we obtain, 

^c(l,'-'^)=[ -Jcm<l,i..)) -(^.„(q,ia,.>\ I 

. X / JsidM^^^u)) Js[l + Js(^]M,^^ 
Ds = , M , V / II . / x\i 72/ II / • X 

\ Js[i + Js (<^c/(q> ^^i^) )J Js {fffi^^ ^^^) 




and equivalent expression for the transverse spin part ^(q, icj;^). The different bub- 
bles {ip{c[,iuju)) in the latter expressions represent the susceptibilities ip{(\,iujy) (cf. ap- 
pendix B) averaged over the distribution of fields cxp[— AS'gjj(r, A, ^j, ^^)]. At that 
point, the problem is highly complex and cannot be solved exactly. Therefore, we in- 
troduce a local and instantaneous approximation which is equivalent to averaging over 
gaussian distributions of local fields (n, Aj, ^j., ^^i) with variances (T\, af and ctc 
defined as 
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a 



J = 3^[D|(q,z^,y, (10) 



and similar expressions for the three other variances. Eq.(p!0[) is obtained from Eq.(|^) by 
letting the summation over (g, iujy) act on the matrix elements only, hence defining the 
local and instantaneous fluctuations of the various degrees of freedom. One expects that 
at low temperatures, the fluctuations of rj, Aj and only bring corrections to their 
saddle-point values. Those corrections are negligible compared to the T-dependence 
brought by the fluctuations of and will not be considered here. So we will neglect in 
the following the fluctuations of r^, Aj and and focus on the fluctuations of only. 
Then {ip{c{,iujy)) is simply given by 

I ( ■ ^^ (/^^0,</'(q'^^^>'^0,Ao,0.>^co)exp[-e|y(2(T2)]) 

= (m,.M-e,j{2.))\) ■ 

Eqs. ( pi][|TT| ) are the basic equations of the paper from which the dynamical spin 
susceptibility is extracted and the critical phenomena around the quantum critical point 
(QCP) are studied. Let us note that, when the averages acting on the various bubbles 
(p{(\,iujy) in Eq.(|D are omitted, one recovers the standard results of the random phase 



approximation. As pointed out by Lonzarich and Taillefer the corresponding self- 
consistent renormalized spin fluctuation SCR-SF procedure can be analyzed within the 
Ginzburg-Landau formalism in which the local free energy is expanded in terms of a 
small and slowly varying varying order parameter m(r) as /(r) = fo + \a |m(r)|^ + 

^6 |m(r)|^ + \cY^ |Vm*(r)|^ + 0(6) and the following approximation is made |m(r)|^ = 

i 

(|m(r)|^) |m(r)|^. Hence (|m(r)|^) is evaluated via the fluctuation-dissipation theorem. 



III. DYNAMICAL SPIN SUSCEPTIBILITY 




Let us now derive the expression of the dynamical spin susceptibility. For that 
purpose, we study the linear response Mf to the source-term — 2S/.B (we consider B 
colinear to the z -axis). The effect on the partition function expressed in Equation (^ ) 
is to change the Hj{t) to H'f{T) 

Hj (r) = Ew [ j j^^^rl, E (-^5^C" - 

\ oi=x,y,z 

+ JcE^^*^^ + JsEi^fAc.. (12) 

Introducing the change of variables = — iB/Js, we connect the /- magnetization 
and the // - dynamical spin susceptibility to the Hubbard Stratonovich field ^j. 

M^__^dLnZ_ ,. 

^ff = ~p dB^dBP ^ ~ ) ^ ^^^^ ) ) ■ ^^^^ 
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Using the expression (||) for the boson propagator D5^(q), we get for the Matsubara 
spin susceptibihty 

X„(^i,^-,^- , , . , <^^(fa7-» (14) 

The diagrammatic representation of Equation ( p!^ is represented in Fig. 2 when the 
averaging is omitted. In the low frequency limit, one can easily check that the dynam- 
ical spin susceptibility may be written in terms of intra- and inter-band suceptibilities 
corresponding respectively to particle-hole pair excitations within the lower a band and 
from the lower a to the upper j3 band. 



X//(q, tu,) = — , . \\ - (15) 

^ ^5 \Aaa iq, ^^^j; \Xap[(l, 



The expressions of the susceptibilities ^^^(q, icu,^) and Xaf^i^^'^^u) corresponding 
to intra- and inter-band particle-hole excitations are given in Appendix B. The latter 



expression is reminiscent of the behaviour proposed by Bernhoeft and Lonzarich |3T|] to 
explain the neutron scattering observed in U Pt^ with the existence of both a " slow" and 
a "fast" component in x (q, uj)/uj due to spin-orbit coupling. Also in a phenomenological 
way, the same type of feature has been suggested in the duality model developed by 
Kuramoto and Miyake [Q. To our knowledge, the proposed approach provides the first 
microscopic derivation from the Kondo lattice model of such a behaviour. Expanding 
the various susceptibilities ip{q,iUu,ro, Xo,^f-,^co) up to the second order in and 
making use of Eq. (pUf ), one can draw from Eq. ([T5|) the following expression of the 
dynamical spin susceptibility (taking the analytical continuation iu!,y ^ uj + i6) 

^^^^^'"^ - l-J|x.«(q,-)x.,(q,-) + Aaf ^^^^ 

where A is a constant which can be evaluated from the expansion of {xaaisiii'^v)) and 
i^^p{(\,iLOy)^ up to the second order in a^. The variance ctj defined by Eq. (p!0D can as 
well be expressed as a function of X//(q) ^^u) 



3 /•+00 

J2 Xffi^^i^u) = -J2 coth—Xjf{ci,uj)duj. (17) 



3 f+oo 

,2 



The two last equations Eqs. ( [Lq - |17|) provide a self-consistent determination of aj and 
hence of the spin susceptibility. We will successively present now (i) in section IV, the 
result for the dynamical spin susceptibility in the random phase approximation (RPA) 
giving information on the nature of the collective mode that is involved; (ii) then in 
section V, the implications of that mode in the critical phenomena around the quantum 
critical point. 
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IV. RPA-DYNAMICAL SPIN SUSCEPTIBILITY NEAR THE 
ANTIFERROMAGNETIC WAVE VECTOR : DISCUSSION ON THE NATURE 

OF THE COLLECTIVE MODE 



As pointed out at the end of section II, the random phase approximation (RPA) 
corresponds to the absence of any averaging in the different equations. Thereby the 
expression of the dynamical spin susceptibility at the RPA level is given by Eq. ([TBI) 
taking a'j = 0. We now discuss the q- and u- dependence of the RPA dynamical spin 
susceptibility around the antiferromagnetic wave- vector Q. 

The bare intraband susceptibility Xqq(Q + ^t, t^) is well approximated at |q| ^ |Q| 
by a lorentzian Xao(Q + c[, oj) = Paa/i^ ~^ i^/'^o + hq^) where paa = Xaa 

(Q,0) and 

Fq is the relaxation rate of order = Tk- As expected for the antiferromagnetic 
case, the relaxation rate remains finite when q goes to zero reflecting the fact that the 
fluctuations of the order parameter are not conserved. The bare interband susceptibility 
Xa^(Q + q, t^) can be considered as purely real and frequency independent equal to pafi 
in the low frequency limit. Fig. 3 gives the continuum of the intraband and interband 
particle-hole excitations x'La ^^^d Xap 7^ 0- -Due to the presence of the hybridization 
gap in the density of states, the latter continuum displays a gap ranging from 2(To, the 
value of the direct gap for the wave- vector 0, to 2 \yF\, the value of the indirect gap at 
Q. The RPA dynamical spin susceptibility is given by 



with 



F(g) = Fo(l-/,) 

where Ig = I — bq^ and / = J^PaaPap (of order 1 near the antiferromagnetic instability). 

The bulk staggered susceptibility X//(Q) diverges at / = 1. The frequency depen- 
dence of X//"* "(Q + Q)"^) is a lorentzian with a vanishing relaxation rate F(0) at the 
antiferromagnetic transition. The corresponding excitation can be analyzed as an an- 



tiferromagnetic paramagnon mode that softens at the magnetic transition. Eq. ([18 
provides us with the dispersion of that mode. The corresponding dynamical exponent z 
is found to be equal to 2 due to the overdamping of the mode when enters the continuum 
of the a-a particle-hole excitations. We will show how the dynamic exponent z strongly 
affects the static critical behavior. This is due to the fact that, contrary to the T > 
phase transitions, statics and dynamics are intimately linked in quantum (T = 0) phase 
transitions. We turn now to the resulting quantum-classical crossover occuring at finite 
temperature. 
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V. CRITICAL PHENOMENA AROUND THE ANTIFERROMAGNETIC 
QUANTUM CRITICAL POINT (AF-QCP) 



The parameter which controls the temperature dependence of the thermodynamical 
variables near the AF-QCP is aj which expresses the thermal local fluctuations of the 



staggered magnetization. Using Eqs. (|T^ and (|r 



aj = 5: 5(Q + q) ~ / 5(Q + q)g^-^dg, (19) 
where ^(Q + q) is the static form factor 

S{Q + q) = - coth(/5^/2)x//Q + q, uj)diu. (20) 

Depending on the temperature-scale considered, quantum or classical behavior is 
observed. In general, the thermal fluctuations aj- of the magnetization should depend 
on the total self-consistent renormalized dynamical spin susceptibility X//"(Q + Q)'^) 
as defined in Eq. ([T5|) . However, in practice, we will content with its truncated form 
Xff^ "(Q + q, cj) derived at the RPA level. The temperature dependence of the static 
form factor ^(Q + q) is given by 

fiT T 
S{Q + q) = —x'ffiQ + q) tan-^(^). (21) 

A first cross-over temperature T/ = 2ro(l — /) appears which separates the quantum 
from the classical regime 

(i) T < Tj : quantum regime. One can show that the static form factor 5'(Q + q) 
exhibits a quadratic temperature dependence 

s(Q + q) = ^^^^^^^ (22) 

TT 1 o(l - iq) 

The thermal local fluctuations cx/^ of the magnetization are also quadratic as a function 
of the temperature. They can be easily evaluated from Eq.(^T]). Taking advantage of 



the fact that the integrand S'(Q + q)g'^ ^ is peaked at gi = y (1 — /)/6, the temperature 
dependence of o"/^ is also found to be quadratic 



The physics is quantum in the sense that the fluctuations on the scale of the magnetic 
correlation length ^ have energy much greater that the temperature ksT. 

(ii) T > Tj : classical regime. In this regime, the static form factor ^(Q + q) shows 
a linear temperature dependence. As long as q is smaller than the thermal cut-off q* 
defined by uJmaxiQ*) = T, S{Q, + q) takes the following form 
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S{Q + q) = 3x'ff{Q + q)T. 



(24) 



The thermal local fluctuations crj^ of the magnetization are deduced from Eq. (pT]). 
The main contribution arises from the integration over [gij^*]. The result depends on 
the dimensionality d. One finds o"/^ ~ T^/^ at = 3 and ctj ~ Tin T at = 2. This 
classical regime corresponds to the case where the energy of the mode on the scale of ^ 
becomes less than ksT. 

Table I recapitulates the expressions that we get for the thermal local fluctuations 
aj of the staggered magnetization depending on the temperature scale and on the di- 
mensionality. Following Eq. (p!6|), the staggered spin susceptibility is given by 

where = Paa+Pap- Incorporating the values of the thermal local fluctuations cx/^ 

reported in Table I, we are able to identify a second cross-over temperature Tu . Above 
Til 1 CT/^ becomes larger than (1 — /) and all the physical quantities are controlled 
by the temperature only. Tn ~ (1 — J)^/^ at c/ = 3 and Tn ~ (1 — /) at (i = 2. 
As expected in the two-dimensional case, the T = transition is at its upper critical 
dimension since z + d = A 'm the antiferromagnetic case characterized by a dynamic 
exponent z = 2. That corresponds to the marginal case for which the intermediate 
regime Ti < T < Tn is squeezed out of existence. The results for the d = 3 and d = 2 
cases are respectively summarized in Figures 4 and 5 which picture the different regimes 
of behaviors reached depending on the values of the temperature T and of the control 
parameter / = J'^paaPap- When / > 1, long-range antiferromagnetic order occurs when 
T is smaller than the Neel temperature T/y. T/v ~ (/ — 1)^^^ at c? = 3 and T^r ~ (/ — 1) 
aid = 2. 



VI. DISCUSSION 

The low temperature physics is dominated by the presence of an antiferromagnetic 
paramagnon mode that softens at the magnetic transition. The results depend crucially 
on the value oi d + z where d is the spatial dimension and z is the dynamic exponent 
associated to that mode. Due to the overdamping of the mode when it enters the 
continuum of intraband particle-hole excitations, the dynamic exponent z is found to 
be equal to 2 in the antiferromagnetic case. As far as the dimensionality is concerned, 
there is no doubt about the dimensionality associated with the heavy quasiparticles in 
those systems that is clearly of 3. However the question about the dimensionality of the 
critical magnetic fluctuations which the quasiparticles are coupled to, is still an open 
problem. Up to now, it has been assumed that the critical fluctuations of the QPT are 
dominated by the three-dimensional antiferromagnetic correlations. This would impose 



a description by a quantum critical theory with d = 3, z = 2. Recently Rosch et al [pq] 
proposed, on the basis of the neutron scattering data in CeCwg-x^^x, that the critical 
magnetic fluctuations are two-dimensional which leads to a QCP with d = 2, z = 2. Let 
us summarize the different regimes of behaviors that we get in each of those two cases 
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{z = 2 , d = 3 OT 2) depending on the values of the temperature T and of the control 
parameter I = J'^PaaPa/s- As one can see, the results that we obtain are very similar to 
those established by Hertz |^ and Millis [|1^ using renormalization group approaches 



in the spin-fluctuation theory. 

{i) Case d = 3. A long-range antiferromagnetic phase occurs when J > 1 below 
the Neel temperature ~ (/ — 1)^/^. A first cross-over temperature T/ ~ (1 — /) 
separates the quantum from the classical regime. In the regime /(/<!) and T < T/), 
the physics is quantum in the sense that the energy of the relevant mode on the scale of 
^ is greater than the temperature ksT. Since d + z > 4, the T = phase transition is 
above its upper critical dimension and the various physical quantities depend upon the 
parameter /. The magnetic correlation length diverges at the magnetic transition ^ ~ 



1/y (1 — J) and the staggered spin susceptibility behaves as x'q = ~-^ + c^^^)• 

Regime II and /// are both classical regimes characterized by large thermal effects since 
the fluctuations on the scale of ^ have energy much smaller than ksT. In Regime II 
(Ti <T < Tjj with Tjj ~ (1 - 1)^/^), the magnetic correlation length ^ ~ - I) 

is still controlled by (1 — /) even though modes at the scale of ^ have energies less 
than ksT. On the contrary, the staggered spin susceptibility is sensitive to the thermal 
fluctuations Xq = - I + aT^'^). In Regime /// (T > Tn), the thermal 

dependence of physical quantities becomes universal and both ^ and Xg are controlled 
by the temperature : ^ ~ and Xq = x^ff{Q)/T^^'^- 

(ii) Case d = 2. Then, since z = 2, z + d = 4: and the T = phase transition 
is at its upper critical dimension. The physics is qualitatively similar to the case d = 
3 with stronger fluctuation effects particularly in the classical regime. A long-range 
antiferromagnetic phase occurs when / > 1 below the Neel temperature T/v ~ (/ — 1). 
For d = 2, the two cross-over temperatures T/ and T// coincide so that Regime II 
of Figure 4 is squeezed out of existence. Regime I {T < Tj with T/ ~ (1 — /)) is 
the quantum regime in which the thermal fluctuations are negligible : ^ ~ 1/ — /) 

and x'q = — I + aT'^). Regime /// (T > Tj) is the unique classical regime 

characterized by very strong fluctuation effects since /c^T is larger than the energy of 
the relevant mode on the scale of ^ . Both the magnetic correlation length and the 



staggered spin susceptibility are then controlled by the temperature : ^ ~ 1/vTlnT 
and XQ = xUQ)/iT\nT). 



VII. CONCLUDING REMARKS 



We considered the 5 = 1/2 Kondo lattice model in a self-consistent one-loop approx- 
imation starting from a generalized Hubbard-Stratonovich transformation of the Kondo 
interaction term. The model exhibits a zero temperature magnetic phase transition at 
a critical value of the Kondo coupling. The transition is usually antiferromagnetic but 
it may be incommensurate depending on the bandstructure considered. The low tem- 
perature physics is controlled by a collective mode that softens at the antiferromagnetic 
transition with a dynamic exponent z equal to 2. 
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A quantum-classical crossover occurs at a temperature Tj related to the charac- 
teristic energy-scale of that mode. Heavy-fermion systems are usually believed to be 
3- dimensional. However since some recent inelastic neutron scattering experiments per- 
formed in CeCuQ-x^'iJ'x show that the critical magnetic fluctuations which the quasipar- 
ticles are coupled to, are 2-dimensional, we both considered the cases d = 2 and d = 3 
with z = 2. The low-temperature behavior of the system is deduced, with predictions 
for the temperature dependence of the physical quantities such as the magnetic corre- 
lation length ^ and the staggered susceptibility x'q- A second crossover temperature 
T// appears above which the thermal dependence of ^ and x'q becomes universal and is 
uniquely controlled by the temperature. It would be very interesting now to study the 
temperature behavior of other physical quantities as the specific heat or the transport 
properties. 

A number of significant issues remain to be addressed. Among others, we will men- 
tion the possibility of getting a different value of the dynamic exponent for instance 
z = 1 if nesting effects are considered corresponding to an absence of damping of the 
mode when located in the gap of excitations around the antiferromagnetic vector. Dis- 
order is expected to play a crucial role in this problem. In that direction, some recent 
works P^ , p!7| pointed out the formation of some "hot lines", i.e., points at the Fermi 
surface linked by the magnetic Q vector. On the hot lines, the quasiparticle scattering 
rate is linear in temperature while, away from the hot lines, it acquires the standard 
Fermi liquid form in T^. At low temperature, the "cold" regions are shown to short- 
circuit the scattering near those lines, and induce a behavior of the resistivity. The 
disorder is found to emphasize the contribution of the hot lines, eventually leading to a 
non- Fermi liquid behavior. It would be worthy to study the effects of the "hot lines" in 
presence of disorder within the critical phenomena description presented in this paper. 

We would like to thank Nick Bernhoeft, Piers Coleman, Jacques Flouquet, Patrick 
A. Lee and Gilbert Lonzarich for very useful discussions on this work. 

* Also at the Centre National de la Recherche Scientifique (CNRS) 
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A. Appendix A: Saddle-Point 



The saddle-point solution is obtained from Eq. (^) for space and time independent 
fields $0; -^0; C/o ^co- III the magnetically-disordered regime {^f^ = C,co = 0), it leads 
to renormalized bands a and (3 as schematized in Figure 1. Noting ctq = Jc^o 
Ef = Ao, aL|0 > ^iid /?fco-|0 > the eigenstates of 



Go-^<^(k,r) 



dr + £:fc cr* 

(To dr+Ef 



(26) 



with respectively the eigenenergies (j)r + and + E'^^ . In the notations: Xk 



£k - Sf, =Ef^ - Ef and Aa 



+ 4(To, we get 



l/± = {xk ± A,) /2. 



(27) 



Let us note Ul^ the matrix transforming the initial basis {c\^ fl^) to the eigenbasis (a^^ 
The Hamiltonian being hermitian, the matrix Uka is unitary : UkaUl^ = UljJka = 

( -Vk Uk ^ 



1. In the notation U, 



t 

ka 



Uk Vk 



we have 



Uk 



k? 



Xk 

1 + -^ 
Afc 



(28) 



Vk 



1 + (^o/i/fc y 



Xk 

Afc 



The saddle-point equations together with the conservation of the number of conduc- 
tion electrons are written as 



iV r 



S ka 



(29) 



^ ka 



TtY^^I npiEk ). 



^ ka 



Their resolution at zero temperature and for ric close to 1 leads to 

lypl = Dexp [-2/ (poJc)] 



(30) 
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/i = 0, 

where yp = fi — ep and po is the bare density of states of conduction electrons (po = 
1/2D for a flat band). Noting y = E — ep, the density of states at the energy E is 
p (E) = Po (1 + (jQ/y^). If Tie < 1, the chemical potential is located just below the upper 
edge of the a -band. The system is metallic. The density of states at the Fermi level 
is strongly enhanced compared with the bare density of states of conduction electrons : 
p{Ep)/pQ = (1 + o'o/y'p) ~ l/(2po IufD- This corresponds to the flat part of the a-band 
in Figure 1. It is associated to the formation of a Kondo or Abrikosov-Suhl resonance 
pinned at the Fermi level resulting of the Kondo effect. The low-lying excitations are 
quasiparticles of large effective mass m* as observed in heavy-Fermion systems. Also 
note the presence of a hybridization gap between the a and the P bands. The direct 
gap of value 2cro is much larger than the indirect gap equal to 2\yp\. The saddle-point 
solution transposes to N=2 the large-N results obtained within the slave-boson mean- 
field theories (SBMFT). 

B. Appendix B : Expressions of the different bubbles 

The expressions of the different bubbles appearing in the expression of the boson 
propagators (cf. Eq.|^) are given here (with i=l, 2, m or ff) 

^j(q, iuJu) = ipi{q, iu^u) + V'i(-q, -i^^u) (31) 
ipi{ci,iuj^) = -'^ G'^^,j(k + q,zcj„ + zu;^)G'^^^,(k,zcj„) 



^ ka.iujn 



1 



kaAuJn 



P I. 



^jM^i^u) = -'^ Y G]^^{k + q,iUn + iuJu)G^fco{Ki^^n), 
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where G^^^(k, iun), Gjjy(k, iun) and Cj^^ik, iun) are the Green's functions at the saddle- 
point level obtained by inversing the matrix G'Q(k, r) defined in Equation (^). 

The different bubbles ipff{q,iuj^), ipcc{(l,iuJu) and ipfc{q,,iuju) can be expressed as a 
function of the Green's functions associated with the eigenoperators al^ and 

Gff{k, iUn) = ulGaaik, iuj^) + WfcG'^/3(k, iujn) (32) 
G^{k, iUn) = vlGc,a{k, iujn) + ulGfSfsik, iUn) 

Gcf{k,iujn) = Gfc{k,iujn) = -UkVk[Gaa{k,iuJn) " G/j/j (k, , 

where Gaa{k, iujn) and G^jjik, iuJn) are the Green's functions associated to the eigen- 
states Ofco-lO >and /^^g-lO >. In the low frequency limit, one can easily check that the 
dynamical spin susceptibility may be written as 

for both the longitudinal and the transverse parts. 



_ rp + 



— / • ^ V^/ 2 2 I 2 2 \ ^ V K / ^ V"/c 

X„^(q, tuj^) = - 2^{UkVi,^g + ffcMfc+gj— — p — 
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TABLE CAPTIONS 



Table I : Predictions for thermal local fluctuations Oj^ of the staggered magnetization 
depending on the temperature T and on the dimensionality d. Tj {\ — 1^ represents 
the first cross-over temperature separating the quantum from the classical regime in the 
vicinity of the antiferromagnetic-quantum critical point (AF-QCP) . 





quantum regime T < T/ 


classical regime T > Tj 


= 3 


a) ~ xW(Q)^Vro 




d = 2 


a) ^ x'ff{Q)TVro 


a} -TlnT 
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kp 

Figure 1: Energy versus wave- vector k for the two bands a. and /3. Note the presence 
of a direct gap of value 2(Jo and of an indirect gap of value 2 
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Figure 2: Diagrammatic representation of Equation ([1^) for the dynamical spin 
susceptibility xjj(q, cu). 
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Figure 3: Continuum of the intra- and interband electron- hole pair excitations 
X^^{q,u!) 7^ and X^^(?, ^ 0. Note the presence of a gap in the interband tran- 
sitions equal to the indirect gap of value 2\yF\ at q — kp, and to the direct gap of value 
2(To at g = 0. 
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Figure 4 : Phase diagram in the plane (T, I) for dimension equal to 3. The shaded 
region represents the long-range antiferromagnetic phase bordered by the Neel temper- 
ature Tjv. The unshaded region marks the magnetically-disordered regimes /, II and 
III associated to different behaviors of the system. Regime / is the quantum regime 
in which the energy of the relevant mode on the scale of ^ is much greater than ksT. 
The magnetic correlation length is ~ 1/ ^J {1 — I) and the staggered spin susceptibil- 
ity is x'q = X^ff{Q)/{^ — I + aT"^). Regime // and /// are both classical regimes in 
which the thermal effects are important since the fluctuations on the scale of ^ have 
energy much smaller than ksT. In Regime //, ^ ~ 1/^(1 — 1) is still controlled by 
(1 — /) but the staggered spin susceptibility is sensitive to the thermal fluctuations : 
Xq = x}fiQ)/i^ — I + aT^^"^). In Regime ///, both ^ and Xq ^-^e controlled by the 
temperature : ^ - l/T^/^ and x'q = X^ff{Q)/T^^^- 
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Figure 5 : Phase diagram in the plane (T, I) for dimension equal to 2. The shaded 
region represents the long-range antiferromagnetic phase bordered by the Neel tem- 
perature T^v- The unshaded region marks the magnetically-disordered regimes I and 
III. In that d = 2 case, the equivalent of Regime II in Figure 3 is squeezed out of 
existence since T/ = T//. Regime / is the quantum regime in which the thermal fluc- 
tuations are negligible : ^ ~ 1/^(1 — /) and x'q = X^ff{Q)/{^ — I + aT'^). Regime /// 
is the unique classical regime in which both the magnetic correlation length and the 
staggered spin susceptibility are controlled by the temperature : ^ ~ l/y/TlnT and 

XQ = x?/(Q)/(rinr). 
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